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''’■The  most  fundamental  environmental  conditions  inhibiting 
laser  use  are  high  rain  rates  and  low  visibility  due  to  fog,  haze 
or  dust.  This  report  calculates  the  frequency  of  occurrence  of 
these  conditions  as  a  function  of  geographic  location.  The 
geographical  locations  chosen  were  those  which  may  have  possible 
laser  operations,  and  include  the  Persian  Gulf,  North  Atlantic, 
Mediterranean,  Southeast  Asia,  Korea,  and  the  Caribbean  Sea. 

With  a  limiting  3  hour  average  rain  rate  of  3  mm/hr,  the 
frequency  of  occurrence  of  the  rain  limitation  is  no  greater  than 
3%.  This  depends  on  latitude  -  higher  latitudes  having  lesser 
rain  rates  -  as  well  as  time  of  year.  The  limiting  visibility 
condition  was  considered  to  be  1  km.  Visibility  less  than  1  km 
was  generally  related  to  fog  and  as  such  is  also  related  to 
latitude.  In  the  tropical  areas  the  visibility  limitation 
occurred  less  than  of  the  time  regardless  of  time  of  year. 

In  polar  and  subpolar  locales,  limited  visibility  occurred  over 
20$  of  the  time  (Kamchatka,  in  July). 
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1.  INTRODUCTION 


1.1  PURPOSE 

The  purpose  of  this  effort  was  to  determine  the  frequency  of 
occurrence  of  weather  conditions  which  preclude  operation  of  the 
high  energy  laser  system  (HELS)  in  real  operational  conditions. 
The  limitations  include  range  obscuration  by  rain,  fog  and  haze, 
as  well  as  the  humidity  and  aerosols. 


1.2  OBJECTIVE 

The  objective  of  this  project  was  to  produce  a  climatology 
of  the  following  variables  for  various  Marsden  Squares:  rainrate; 
visibility;  total  extinction  coefficient  over  the  wavelength  band 
used  by  the  HELS  poi nter-tracker ( 8-1 2  pm);  molecular  extinction 
coefficient  for  the  mid-range  infrared  advanced  chemical  laser 
(MIRACL)  spectrum.  The  climatology  was  expressed  in  terms  of  the 
average  frequency  with  which  specified  values  of  the  above 
variables  were  reported. 


1.3  BACKGROUND 

The  decision  to  deploy  a  High  Energy  I.aser  System  (HELS)  as 
part  of  a  weapon  suite  depends  on  the  vulnerability  of  the  HELS 
to  expected  weather  conditions.  The  environmental  factors  which 
influence  the  system  are  scattering  and  absorption  by  atmospheric 
aerosols  and  molecules,  nonlinear  effects,  such  as  thermal 
blooming,  and  turbulence-induced  laser  beam  wander.  These  effects 
are  not  directly  measurable  with  standard  meteorological  instru¬ 
mentation,  but  correl  ati on  model  s  can  be  used  to  relate  archived 
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meteorological  measurements  to  the  desired  meteorological 
eters  .  The  collected  and  analyzed  results  predict  the  frequency 
of  occurrence  of  weather  conditions  affecting  the  HELS. 


2.  DATA 


2.1  DATA  BASE 

A  consolidated  Data  Set  (CDS)  has  been  created  at  Fleet 
Numerical  Oceanography  Center  (FNOC)  using  the  National  Climatic 
Center's  marine  surface  observation  archives  (TDF-11)  and  FNOC's 
surface  observation  archives  (SPOT).  These  data  have  been 
organized  by  Marsden  Square  and  are  stored  on  magnetic  tape. 

Each  tape  consists  of  multiple  files  of  Marsden  Squares.  Within 
each  Marsden  Square  data  are  ordered  by  1-degree  subsquare. 

Within  each  subsquare  data  are  ordered  chronologically  over  the 
period  1946-1977, 

2.2  SOURCE  DATA 

Data  used  in  the  analysis  was  a  subset  of  the  CDS.  The 
period  1963-1973  inclusive  was  chosen  as  the  analysis  period. 
Therefore,  up  to  eleven  years  of  observations  have  been  used  to 
calculate  a  set  of  environmental  statistics. 

The  analysis  concentrated  on  geographic  areas  which  were 
considered  possible  future  HELS  operating  areas.  These  areas 
were  the  Persian  Gulf  (Marsden  Sq.  102-103),  North  Atlantic 
(Marsden  Sq.  145-147),  Southeast  Asia  (Marsden  Sq.  26),  Korea 
(Marsden  Sq.  131-132),  and  the  Caribbean  (Marsden  Sq.  44).  These 
locations  are  described  in  Figure  1  and  Table  1. 

As  each  record  was  processed  the  data  were  checked  for  (1) 
parity  errors,  (2)  correct  packing  arrangement  of  the  observa¬ 
tions  (i.e.,  the  number  of  words  in  the  record  must  be  equal  to 
some  multiple  of  5),  and  (3)  agreement  of  the  air- sea  temperature 
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Marsden  square  in  which  the  sta1!r'r<  !s  situated  at  the  time  of  observation 


difference  with  the  difference  of  air  temperature  and  sea  surface 
temperature.  Records  which  failed  these  checks  were  dropped  from 
the  analyzed  data  base. 


Table  1.  List  of  Marsden  Squares  analyzed 
Square  Geographic  Area 


102-103 


Persian  Gulf 


AraDian  Sea 
Southeast  Asia 


131  -132 


Korea 


Kamchatka  Peninsula 
Caribbean  Sea 


145-147 


143-144 


North  Atlantic  Ocean 
Norway 

Mediterranean  Sea 


2.3  DISTRIBUTION  OF  OBSERVATIONS 

A  wide  variation  in  the  temporal  distribution  of  the  observa¬ 
tions  was  noted.  A  given  month  of  some  years  had  few  or  no  reports 
while  the  same  month  of  other  years  had  many  reports.  The  temporal 
distribution  varied  for  each  Marsden  Square  analyzed.  It  was 
determined  that  a  simple  averaging  of  the  data  over  the  entire 
analysis  period  would  yield  mean  values  that  would  be  biased  by  the 
data  from  years  with  many  observations.  An  alternate  method  of 
computing  mean  values  was  chosen  in  which  monthly  averages  were 
calculated  for  each  year  with  one  or  more  reports  for  the  given 
month.  These  monthly  averages  were  then  averaged  to  obtain  mean 
values  for  the  entire  analysis  period.  This  method  of  computing 
the  means  had  the  effect  of  limiting  the  bias  introduced  by  years 
that  had  many  reports  relative  to  other  years.  The  averages  of 
the  parameters  were  averaged  over  the  Marsden  subsquare  averages 
to  produce  a  Marsden  square  average.  The  procedure  for  calculating 
mean  values  for  each  Marsden  square  gave  equal  weight  to  all 
subsquares  with  reports. 
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3.  METHODOLOGY 


3.1  OVERVIEW 

The  parameter  correlation  models  described  in  this  section 
were  used  to  calculate  the  HEL's  characteristics  associated  with 
each  observation  recorded  on  the  CDS  tapes.  The  resulting  values 
for  the  HEL's  variables  were  grouped  by  month  and  Marsden  subsquare 
and  averaged.  The  resulting  means  then  were  grouped  by  Marsden 
square  for  the  entire  analysis  period  and  averaged. 

3.2  COMPUTER  PROCESSING 

A  computer  program  was  developed  to  process  the  thousands  of 
weather  observations  that  formed  the  data  base  for  the  analysis. 

The  data  recorded  on  the  CDS  tapes  was  known  to  be  relatively  free 
of  data  gaps  and  errors  so  only  a  few  quality  control  checks  were 
included  in  the  program.  Speed  in  operation  and  the  production 
of  conveniently  formatted  output  (tables  and  graphs)  were  the 
main  features  of  the  program. 

Below  is  a  list  of  the  primary  functions  of  MARPLT  (Marsden 
Square  Plotting)  program: 

1.  Read  the  data  from  the  CDS  tape. 

2.  Establish  that  the  record  is  acceptable  for  processing 
(check  for  parity  error,  EOF  marker,  packing  arrangement  of  the 
observations ) . 

3.  Unpack  the  coded  data  using  a  masking  expression  and  the 
SHIFT  intrinsic  function. 

4.  Select  only  those  observations  taken  during  the  1963- 
1973  time  period. 
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5.  Calculate  values  for  the  following  variables:  visibility, 
rainrate,  molecular  absorption  coefficient,  molecular  extinction 
coefficient,  and  extinction  coefficient  for  the  MIRACL  spectrum. 

6.  Calculate  the  cumulative  frequency  of  occurrence  of: 
visibility  less  than  or  equal  to  1  km,  rainrate  greater  than  or 
equal  to  2  millimeters  per  hour  (mm/hr),  molecular  extinction 
coefficient  (8-12  pm)  greater  than  or  equal  to  0.8  per  kilometer, 
and  MIRACL  spectrum  extinction  coefficient  greater  than  or  equal 
to  0.28  per  kilometer. 

7.  Calculate  the  frequency  of  occurrence  of  rainrates  for 
21  intervals  ranging  from  0  to  greater  than  10  mm/hr  in  0.5  mm/hr 
increments  and  grouped  by  month. 

8.  Calculate  the  cumulative  frequency  of  occurrence  of 
rainrates  for  the  21  intervals  for  the  months  of  March,  June, 
September,  and  December. 

9.  Calculate  the  mean  and  standard  deviation,  grouped  by 
month  and  Marsden  subsquare,  of  the  following  parameters:  molecular 
extinction  coefficient  (8-12  pm),  and  MIRACL  spectrum  extinction 
coef  fi  ci  ent . 

10.  Produce  an  output  of  tables  and  graphs  summarizing  the 
statistics. 

The  flow  chart  in  Figure  2  illustrates  the  flow  of  data  through 
MARPLT,  showing  a  closed  loop  which  will  analyze  all  observations 
within  a  Marsden  square  file  before  terminating.  The  numbers 
appearing  at  the  upper  right  corner  of  each  box  relate  the  opera¬ 
tion  to  one  of  the  ten  functions  outlined  above.  Details  of  the 
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execution  of  each  of  the  ten  primary  functions  is  given  in  the 
following  sections.  The  parenthetical  numbers  following  each 
heading  indicate  the  relationship  to  the  above  functions. 


3.3  DATA  INPUT  AND  VERIFICATION  (1,2) 

A  buffer  input  statement  was  used  to  transfer  data  from  a 
file  to  a  block  of  memory  one  record  at  a  time.  Then  the  status 
of  the  buffer  input  operation  was  checked  for  an  EOF  or  parity 
error  condition.  If  a  parity  error  was  noted  the  record  was 
discarded.  If  an  EOF  was  encountered,  end-product  statistics  for 
the  Marsden  square  were  calculated.  If  no  parity  error  or  EOF 
was  noted,  the  number  of  words  in  the  record  was  determined.  If 
the  number  of  words  was  a  multiple  of  5,  the  packing  arrangement 
of  the  record  was  acceptable.  Otherwise,  the  record  was  discarded 
Acceptable  records  were  processed  observation  by  observation.  The 
above  procedure  was  repeated  until  all  observations  for  a  Marsden 
square  were  analyzed  and  an  EOF  was  encountered. 

3.4  DATA  DECODING  (3) 

Each  observation  was  first  unpacked  from  a  binary  form.  Each 
observation  was  then  checked  for  data  source  (SPOT  or  TDF-11). 

The  appropriate  scaling  for  pressure  was  applied  depending  on  the 
data  source.  Air  temperature,  dew  point  and  sea  surface  tempera¬ 
ture  were  corrected  for  sign.  Latitude  and  longitude  were  decoded 
with  the  convention  that  positive  (negative)  latitude  was  North 
(South).  Positive  (negative)  longitude  was  east  (west).  The 
remaining  data  were  decoded  with  scaling  and  bias  factors  described 
in  Mil  (1  975  )  . 
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3.5  ANALYSIS  PERIOD  (4) 

Only  those  observations  taken  during  the  eleven  year  period 

1963-1973  were  analyzed.  Selection  of  the  desired  observations 
required  a  simple  check  that  the  year  of  the  observation  was  within 
5  years  of  1968. 

3.6  METEOROLOGICAL  CALCULATIONS  (5) 

Visibility  values  were  taken  directly  from  the  CDS  tapes. 

The  four  remaining  variables  were  calculated  using  correlation 
models  and  the  reported  weather  data. 

3.6.1  Rai nrate 

The  rainrete  was  evaluated  using  the  relation  of  the  present 
weather  code  (IPW)  to  rainrate,  originally  derived  by  Tucker  (1971) 
and  subsequently  modified  by  Bourke  and  Dorman  (1977).  These 
derivations  were  specifically  derived  to  determine  three  hour 
precipitation  tota.s,  rather  than  rainrates.  This  work  infers  an 
average  precipitation  rate  from  the  three  hour  total.  It  must  be 
stressed  that  this  rate  is  often  less  than  the  possible  instanta¬ 
neous  rainrate  within  the  three  hour  period  of  observation. 

3.6.2  Molecular  Extinction  Coefficient 

The  molecular  extinction  coefficient  was  calculated  for  the 
thermal  infrared  wavelength  band  (8-12  ym).  The  extinction 
coefficient  was  calculated  using  a  polynomial  fit  to  the  LOWTRAN  3B 
molecular  transmittance  code.  In  translating  from  transmittance 
to  extinction  a  10  km  range  was  used  to  approximate  the  range 
dependence  of  the  extinction  coefficient. 

11 
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3.6.3  MIRACL  Spectrum  Absorption  Coefficient 


Molecular  absorption  for  the  proposed  MIRACL  spectrum  was 
evaluated  using  an  algorithm  developed  by  Science  Applications, 

Inc.  (SAI).  The  absorption  coefficient  for  each  observation  was 
determined  through  the  solution  of  a  polynomial  which  is  a  function 
of  water  vapor  pressure  (Torr)  and  temperature  (K)  and  has  the 
following  form: 

f(T,P)  =  a,  +  a2T  +  a3P  +  a4P2  +  a$TP2 
+  VP  +  a7T2  +  Vt2  +  a  9 p2?2 

where  T  is  temperature  and  P  is  vapor  pressure.  Values  for  the 
coefficients  a-,  i=  1,2,. ..,9  were  determined  by  SAI  using  a  least 
squares  fitting  procedure  and  are  listed  below: 


a1  =  9 .4351 *1 0"2 

a6  = 

-1  .041 7*1  O'  3 

a2  =  -3  .2344*1  0"4 

II 

2 .6804*1 0"7 

a  3  =  1  .  7890*1  0_1 

a8  = 

1  . 5392*1  O'6 

a4  =  1  .0080*1 0" 3 

a9  = 

1  .2311 *10"8 

a5  =  -7  .0357*1  0"6 

The  polynomial  is  valid  for 

temperatures  ranging  from  240  K  to 

320  K  and  water  vapor  pressures  from  0  to  79  Torr. 

3.7  STATISTICAL  CALCULATIONS  (6, 7, 8, 9) 

The  magnitude  of  the  values  calculated  above  was  checked 
against  specified  limits.  If  these  threshold  values  were  met  or 
exceeded  the  appropriate  frequency  of  occurrence  counters  were 
incremented.  Other  counters  recorded  the  number  of  observations 
for  a  given  month  and  the  number  of  observations  for  a  given 
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Marsden  subsquare.  The  computed  values  for  the  two  attenuation 
coefficients  were  summed  before  calculating  the  mean  and  standard 
deviation  of  the  coefficient  data. 

Data  on  the  CDS  tape  are  chronologically  ordered  within  each 
Marsden  subsquare.  The  processing  described  above  continued  until 
all  observations  for  a  given  month,  year  and  subsquare  were 
analyzed.  Then  program  control  was  transferred  to  the  subroutine 
MNAVG.  When  program  control  was  returned  to  the  main  program, 
the  processing  of  observations  for  the  next  month  but  same  year 
ana  subsquare  was  begun.  (If  observations  from  December  were  the 
last  data  analyzed,  the  month  and  year  would  have  changed  when 
program  control  reverted  to  the  mai n *  program . ) 

After  all  observations  for  a  given  Marsden  subsquare  were 
processed,  program  control  was  transferred  to  the  subroutine 
SUBAVG. 

3.7.1  Monthly  Averaging 

The  subroutine  MNAVG  grouped  information  by  month  and 
performed  some  analysis  of  the  data  accumulated  in  the  main 
program.  No  end-product  statistics  were  calculated  in  MNAVG. 
Values  computed  in  MNAVG  were  based  on  observations  from  a  single 
month,  year  and  subsquare. 

Using  visibility  as  an  example,  the  frequency  of  occurrence 
for  a  given  month  and  year  was  calculated  as  follows: 

FVIS1  =  NVIS/NDMN 

where : 

F  V  IS 1  is  the  frequency  of  occurrence  of  visibility 
<  1  km  for  a  given  month. 
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NVIS  is  the  number  of  observations  that  had  a 


reported  visibility  less  than  or  equal  to  1  km, 
NDMN  is  the  total  number  of  observations  reported 
for  the  month. 

The  results  of  the  calculations  were  identified  by  month 
and  grouped  with  the  results  based  on  data  from  the  same  month 
of  other  years.  (Recall  that  all  of  this  data  is  from  the  same 
Marsden  subsquare.)  Newly  computed  values  of  frequency  of 
occurrence  were  summed  with  previously  computed  values  each  time 
MNAVG  was  called.  The  summed  totals  for  th-^  two  attenuation 
coefficients  were  similarly  grouped  by  month  and  summed  over 
the  entire  data  analysis  period. 

3.7.2  Marsden  Subsquare  Averages 

In  the  subroutine  SUBAVG,  the  various  summations  computed 
in  MNAVG  were  averaged  to  calculate  the  means  and  standard  devia¬ 
tions  of  the  attenuation  coefficients  and  the  mean  frequencies  of 
occurrence  described  previously.  These  statistics  were  grouped 
by  month  and  subsquare.  Using  visibility  as  an  example,  the 
process  was  as  follows: 

1  N 

FVIS2  =  i  l  FVIS1. 
i  =1 

where : 


FV  IS2  is  the  mean  frequency  of  occurrence  of  visi 
bility  £  1  km  in  the  subsquare  for  a  given 
month  over  the  entire  analysis  period. 


y  FVIS1.  is  the  summation  of  the  individual  fre- 
i  =1  1 

quencies  of  occurrence  computed  for  a  given 
month  and  year, 

N  is  the  number  of  years  that  had  reports 
for  the  given  month  and  had  a  possible 
range  of  0  to  11. 

These  mean  subsquare  values  were  then  summed  and  eventually  a 
mean  value  for  the  entire  Marsden  Square  was  calculated. 

3.7.3  Marsden  Square  Averages 

When  all  data  for  a  Marsden  Square  had  been  processed,  the 
program  encountered  an  EOF.  Program  control  was  transferred 
to  a  different  section  of  the  main  program  and  the  summations 
of  subsquare  means  were  averaged  to  produce  the  required  mean 
frequencies  of  occurrence  fo  the  entire  Marsden  Square.  Again 
using  visibility  as  an  example,  the  computations  were  done  as 
f ol 1 ows : 

,  N 


FVIS3 


K  I  FVIS2i 

"  i=l  1 


where : 


FVIS3  is  the  mean  frequency  of  occurrence  of  visi¬ 
bility  <  l  km  in  the  Marsden  Square  for  a  given 
month  over  the  entire  analysis  period, 

N 

l  F V  IS 2 ^  is  the  summation  of  the  mean  frequencies  of 
occurrence  of  visibility  <  1  km  in  each  sub¬ 
square  with  reports, 
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N  is  the  number  of  subsquares  that  had  reports 
during  the  analysis  period. 

The  processing  described  in  this  and  the  preceeding  sections 
can  be  summarized  as  follows: 

*  Monthly  means  were  calculated  for  each  E-O-met  variable 

*  The  monthly  means  were  averaged  to  produce  mean  values 
for  the  subsquares 

*  The  subsquare  means  were  averaged  to  produce  mean  values 
that  were  grouped  by  montn  but  applied  to  the  entire 
Marsden  Square  and  analysis  period. 

This  method  of  computing  the  average  of  means  had  che  effect 
of  limiting  the  bias  introduced  by  years  and  subsquares  that 
had  many  observations  relative  to  others. 

3.8  FORMATTED  OUTPUT  (10) 

Output  from  the  computer  program  included  tables  organized 
by  month  which  listed  the  statistics  described  in  previous  sec¬ 
tions.  Table  2  is  a  sample  of  the  tabular  output.  Histograms 

and  graphs  of  cumulative  frequency  for  the  rainrate  data  were 
plotted  as  shown  in  Figure  3. 

The  output  of  MARPLT  contained  the  mean  latitude,  mean 
longitude,  standard  deviation  of  the  latitude  and  standard 
deviation  of  the  longitude  of  the  Marsden  subsquares  analyzed 
for  each  Marsden  Square.  Subsquares  with  no  reports  were  iden¬ 
tified  by  subsquare  number.  The  number  of  subsquares  with  re¬ 
ports  and  the  number  of  observations  analyzed  for  the  Marsden 
Square  during  the  1963-1973  period  were  given.  The  output  also 
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CUMUL 


F i gu'.'e  3.  Frequency  (top)  and  cumulative  frequency  (bottom) 
of  occurrence  of  rainrates  for  March  for  Marsden  Square  145 


listed  the  identifying  number  of  each  Marsden  Square  analyzed 
and  the  quadrant  of  the  globe  in  which  the  square  is  located. 

All  standard  deviations  were  calculated  using  N  weighting. 

The  subroutine  HISTPLT  contained  the  code  used  to  produce 
the  graphs  on  the  Varian  plotter.  Extended  Core  Storage  (ECS) 
was  required  for  the  plotting  operation.  The  graphs  were  pro¬ 
duced  by  a  simple  bin  sort.  A  total  of  21  intervals  of  rainrate 
were  specified,  ranging  from  0  to  greater  than  10  mm/hr  in  0.5 
mm/hr  increments.  Referring  to  Figure  3(a),  the  histogram 
should  be  interpreted  as  the  frequency  of  occurrence  (percent) 
of  rainrates  falling  into  each  interval.  Reports  were  included 
in  a  bin  when  the  observed  rainrate  was  greater  than  or  equal 
to  the  lower  bin  boundary  value  and  less  than  the  higher  bin 
boundary.  All  observations  of  rainrate  greater  than  or  equal 
to  10  mm/hr  were  grouped  together  and  displayed  in  the  21st 
interval.  For  this  interval  the  histogram  should  be  interpreted 
as  the  frequency  of  occurrence  (percent)  of  rainrates  greater 
than  or  equal  to  10  mm/hr. 

A  variable  vertical  scale  was  used  to  provide  maximum  resolu¬ 
tion  of  the  displayed  data.  The  three  vertical  scales  incorporated 
into  the  program  had  maximums  of  5,  10,  and  20  percent.  Data 
for  the  interval  0.0  to  0.5  mm/hr  were  not  plotted  since  the  over¬ 
whelming  majority  of  observations  fell  into  this  category.  Any 
attempt  to  plot  on  the  same  graph  the  data  for  all  of  the  inter¬ 
vals  would  have  lead  to  scaling  problems  and  difficulty  in  using 
the  plots.  Values  for  the  data  associated  with  the  first  inter¬ 
val  were  listed  in  the  tabular  output. 
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The  cumulative  frequency  of  occurrence  (percent)  corresponding 
to  a  particular  bin  boundary  value  was  derived  by  summing  all  bin 
elements  less  then  the  bin  boundary  value  and  dividing  this  sum  by 
the  total  number  of  bin  elements.  Referring  to  Figure  3(b),  the 
graph  should  be  interpreted  as  the  frequency  (percent)  with  which 
a  rainrate  less  than  the  indicated  value  was  reported.  Again,  one 
of  three  vertical  scales  was  used  to  enhance  the  useability  of  the 
graph  . 


4.  RESULTS 

The  collected  statistics  of  rainfall,  detector  total  extinc¬ 
tion  coefficient  and  MIRACL  molecular  extinction  are  included  in 
Appendix  A.  These  are  in  the  form  described  in  Section  3.8. 

Within  each  geographical  area  statistics  are  provided  for  each 
month.  In  addition  frequency  statistics  of  visibility  and  rain- 
rate  exceeding  their  critical  values  are  provided  for  each  month. 

Program  Mstings  of  programs  relating  archived  meteorological 
data  to  rainrates  and  extinction  coefficients  are  provided  in 
Appendix  B. 


5.  CONCLUSIONS 


The  environmental  conditions  limiting  use  of  a  high  energy 
laser  system  are  basically  high  rain  rate  and  low  visibility. 

These  conditions  limit  the  use  of  the  HEL  system  in  all  regions  ex¬ 
cept  Korea  less  than  6%  of  the  time.  In  the  Korean  area,  visi¬ 
bility  is  limited  by  frequent  fog  conditions.  These  low  visi¬ 
bility  conditions  are  less  frequent  in  winter  (below  2%)  than 
in  summer  when  they  occur  about  10%  of  the  time. 

The  highest  average  rain  rate  limitation  occurs  in  the  tropi¬ 
cal  zones  (for  example  S.E.  Asia  when  rain  limits  laser  use  5% 
of  the  time  in  November).  We  again  stress  the  limitations  cf 
the  rain  rate  algorithm  however.  The  algorithm  is  accurate  in 
northern  latitudes  which  have  long  periods  of  relatively  light 
rain.  The  tropics,  however,  often  have  short  periods  of  very 
intense  rain,  exceeding  25  mm/hr.  The  averaging  techniques  used 
here  "smears"  such  intense  occurrences  over  an  entire  three 
hour  period.  An  average  rain  of  3  mm/hr  in  three  hours,  may 
just  as  well  mean  that  the  entire  rain  fell  in  a  half  hour,  and 
is  characterized  by  a  rain  rate  of  18  mm/hr  for  a  short  period. 
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APPENDIX  A 


TABLES  AND  GRAPHS 


The  data  in  this  appendix  are  presented  in  two  groups:  first, 
tables  of  rainrates  and  bad  weather  statistics;  and  second,  paired 
graphs  of  rainrate  frequencies  of  occurrence  (top)  and  cumulative 
frequencies  of  occurrence  (bottom).  0..a  table  and  eight  graphs  are 
included  for  each  of  the  14  Marsden  Squares  analyzed. 

The  graphs  summarize  their  respective  rainrate  statistics  for 
the  months  of  March,  June,  September,  and  December,  in  that  order. 
These  months  were  chosen  to  represent  the  four  climatic  seasons  of 
spring,  summer,  fall,  and  winter. 

Within  each  of  the  two  groups,  t ~ b 1 es/graphs  are  ordered  by 
Marsden  Square  number  in  the  sequence  shown  below. 
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APPENDIX  B 


MARPLT  PROGRAM  LISTINGS 


This  appendix  contains  program  listings  of  programs  relating 
archived  meteorological  data  to  rain  rates  and  extinction  coefficients. 
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